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The low-energy Hamiltonian used to describe the dynamics of the helical Dirac fermions on the 
surface of a topological insulator contains a subdominant non-relativistic (Schrodinger) contribu¬ 
tion. This term can have an important effect on some properties while having no effect on others. 
The Hall plateaus retain the same relativistic quantization as the pure Dirac case. The height of the 
universal interband background conductivity is unaltered, but its onset is changed. However, the 
non-relativistic term leads directly to particle-hole asymmetry. It also splits the interband magneto¬ 
optical lines into doublets. Here, we find that, while the shape of the semiclassical cyclotron reso¬ 
nance line is unaltered, the cyclotron frequency and its optical spectral weight are changed. There 
are significant differences in both of these quantities for a fixed value of chemical potential or fixed 
doping away from charge neutrality depending on whether the Fermi energy lies in the valence or 
conduction band. 

PACS numbers: 78.20.Ls, 78.67.-n, 71.70.Di, 72.80.Vp 

I. INTRODUCTION 

Recently, the metallic states that exist at the sur¬ 
face of a topological insulator (TI) have been extensively 
studiedip— . These display novel physical properties which 
could possibly be exploited in device applications such 
as quantum computing^. This has lead to a flurry of ac¬ 
tivity both experimentally and theoretically which has 
provided new understanding. Many new materials have 
been discovered^iiS; however, there remains a need to in¬ 
crease the bulk band gap so as to more clearly separate 
the metallic surface states from the bulk insulating spec¬ 
trum. The surface states possess an odd number of Dirac 
points. Here, for simplicity, it is sufficient to consider a 
single Dirac cone centred at the T point of the surface 
Brillouin zone as in Bi 2 Se 3 r^. The Hamiltonian which 
governs the dynamics of the surface fermions involves real 
spini^ in contrast to graphene and other related materials 
where instead pseudospin is involvedi^^—. This is due to 
a strong spin-orbit coupling which provides a dominant 
linear-in-momentum Dirac contribution to the Hamilto¬ 
nian with Fermi velocity vp- There is also a smaller 
non-relativistic contributioni^, which is quadratic in mo¬ 
mentum and characterized by a Schrodinger mass m. In 
this paper, we are particularly interested in understand¬ 
ing the changes in the magneto-optical properties of the 
helical surface fermions that result from the introduction 
of the mass termiS. 

The Schrodinger mass term is a prominent feature of 
the low-energy continuum k ■ p Hamiltonians^id^ which 
form the basis of our understanding of the TI surface 
band structure. Indeed, such work was critical in the dis¬ 
covery of three-dimensional (3D) TIs. These simplified 
Hamiltonians are based on density-functional electronic- 
structure calculations, as first discussed by Zhang et al^ 
and elaborated upon in Refi^^. For the specific cases of 


Bi 2 Te 3 and Bi 2 Se 3 , the band masses are m = 0.16me 
and O.OQrUe, respectively, where rUe is the bare electron 
mass. Thus, for a magnetic field of I Tesla, the respec¬ 
tive Schrodinger magnetic energy scales are ~ 1.25 meV 
and 0.7 meV—. The corresponding Dirac energy scale 
(due to the linear-in-momentum relativistic part of the 
Hamiltonian) is of order 10 meV (a factor of ten larger 
than the non-relativistic contribution). The two scales 
become comparable when the magnetic field B is 40T. 
This results from the Schrodinger terms linear-in-H de¬ 
pendence compared to the \/B scaling of the relativistic 
piece. Therefore, it is important to retain the quadratic- 
in-momentum part of the Hamiltonian when describing 
the surface states. It has been included in the extensive 
literatur o^d°i^°~ — (both theoretical and experimental) 
that exists on the phase offset seen in quantum oscilla¬ 
tions of the Shubnikov-de Haas and de Haas-van Alphen 
effects. There are numerous other optical studies in TIs 
including many in the infrared^^^— ; the most relevant 
for this work are those in the THz rang o^^i^° . Here, 
we mention specifically that of Wu et al^ which pro¬ 
vides information on the cyclotron frequency itself and, 
consequently, directly impacts the present work. After 
the subtraction of a phonon and a lattice polarization 
term, the cyclotron resonance peak is extracted. It was 
found that the quadratic-in-momentum term is essential 
in obtaining a consistent fit to this cyclotron resonance 
data^i. In another study of cyclotron resonance data, 
the Schrodinger pieced was found to be dominant in 
strained films of HgTe (an intrinsic TI with negligible 
bulk carriers^). 

Initially studied TIs where not truly insulating in the 
bulk and the topologically protected surface states could, 
thus, be affected by undesirable mixing with the bulk 
charge carriers. Many strategies have been used to sep¬ 
arate out the surface fermions. In the magneto-optical 
study of Bio. 9 iSbo.o 9 by Schafgans et alS^, several Lan- 
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dau level (LL) series are observed. They all display en¬ 
ergy dependences that vary as the square root of the 
applied magnetic field and none are found to originate 
from the bulk bands since they do not extrapolate to the 
band gap energy. In Bii_ 2 ;Sba;, the observed quantum 
oscillations correspond to a two-dimensional Fermi sur¬ 
face and hence, are identified with the surface states^ii^. 
They also observe a contribution form a 3D Fermi sur¬ 
face which they assign to coupling between the bulk and 
surface states. 

It has also been possible to separate the surface from 
the bulk using time-resolved techniques. The surface 
and bulk states’ transient THz optical response has 
been found to correspond to very different time scales^^. 
Crepaldi et alM- used time-resolved angular-resolved 
photoemission spectroscopy (ARPES) to differentiate be¬ 
tween the two types of states. Such techniques have re¬ 
cently been used^ to show that the Aig phonon mode in 
Bi 2 Se 3 is coupled to both the bulk and surface fermions 
and that this phonon softens on the surface. There are 
many works on the electron-phonon interaction in TIs 
(see, for example, Refsi^^— ). High resolution ARPES 
is a particularly useful technique for this kind of inves¬ 
tigation. Kinks in the dispersion curves of the charge 
carriers are interpreted as coupling to possible bosons. 
Some experiments^ have given large values for the cou¬ 
pling while others^S give much smaller values. In one 
study, the coupling was found to be even smaller and in¬ 
terpreted to imply the topological surface states are not 
only protected from elastic scattering off impurities but 
also from scattering off low-energy phonons^. Such ef¬ 
fects are beyond the scope of this paper. 

Work has also been done to entirely suppress the metal¬ 
lic character of the bulk states. Eor example, this can 
be achieved by doping Cu in Bi 2 Se 3 ^>^. One can then 
study surface fermions without bulk influences. This is 
also the case in Tii_a,Bii+a,Se 2 as documented in the 
work of Eguchi et al^ and Kuroda et al^. 

In the pure relativistic case, there is particle-hole sym¬ 
metry between the conduction and valence bands which 
are both described by the same perfect Dirac cone. The 
introduction of a Schrodinger mass distorts these cones. 
The conduction band narrows with increasing energy 
while the valence band fans out as the energy is lowered. 
This yields an hourglass shape to the overall band struc¬ 
ture. The band bending leads to particle-hole asymme¬ 
try. In this case, doping with holes rather than electrons 
will be different than the pure relativistic case. Here, we 
will emphasize these differences. 

When a magnetic field is applied perpendicular to the 
surface of the TI, Landau levels form. In the pure rela¬ 
tivistic system, their energies are characterized by the 
single Dirac magnetic energy scale Ei = hvF\/eB/h 
where e is the elementary charge. When we account for 
the subdominant non-relativistic contribution, a second 
Schrodinger magnetic energy scale Eq = heB/m enters 
and modifies the eigenenergies of the LLs. In particu¬ 
lar, the level at zero energy which is so important in 


the relativistic case is moved to positive energy {Eq/2). 
While for small values of H, Ex is always larger than 
Eq. These two scales become comparable in size as B 
is increased. This leads to several known effects such as 
the splitting of the interband absorption lines in the lon¬ 
gitudinal conductivityi^i^i. At the same time, the Hall 
plateaus in the transverse DC conductivity retai n^^i^^d^ 
their relativistic quantization which is distinct from the 
non-relativistic values which arise when heB/m is the 
dominant scaled. 

Even when B is small and the semiclassical regime is 
considered, for which the chemical potential |/i| is much 
greater than both Ei and Eq, important modifications to 
the cyclotron resonance frequency (wcr) and its optical 
spectral weight (Wcr) arise when Eq is introduced. Here 
we provide a simple analytic formula for the resulting line 
shape and study the particle-hole asymmetry in both Wcr 
and WcR- While it is the intraband optical transitions 
which are central to the semiclassical regime, the evo¬ 
lution of the interband transitions in the same limit of 
1^1 El and Eq are also of interest. The magnitude of 
the universal background which emerges is completely in¬ 
dependent of the Schrodinger mass but its onset is nol4i. 
In particular, in this work we emphasize the particle-hole 
asymmetry of the onset. 

The paper is organized as follows: In Sec. H, we in¬ 
troduce the Kubo formula for both the longitudinal and 
transverse (Hall) magneto-optical conduct ivityiSii^. We 
also present results for both Recra;a;(D) and Re(Ja;y(D) 
at several values of B and feature their evolution as B 
is decreased sufficiently that the LL structures become 
smooth and reveal a background. In Sec. HI, we start 
from the finite B formulas for the real part of the lon¬ 
gitudinal conductivity and show how the semiclassical 
regime is obtained. The end result is analytic formulae 
for the cyclotron resonance lineshape, frequency and op¬ 
tical spectral weight. We retain the appropriate limit 
as Eq —)> 0. We study the approach to the semiclassical 
limit and discuss the particle-hole asymmetry that arises. 
We compare both a constant |^| for /r ^ 0 and constant 
doping level n for electrons and holes. In Sec. IV, we pro¬ 
vide a discussion of the universal interband background^ 
and its particle-hole asymmetry. An alternative deriva¬ 
tion, based on the finite B formula for the conductivity, 
is provided for the interband background which is addi¬ 
tional to the intraband cyclotron resonance line. Details 
about the asymmetry of the absorption process onsets 
are provided. Finally, in Sec. V, we give a discussion of 
our new results and draw conclusions. 


II. FORMALISM 

The model low-energy Hamiltonian on which this work 
is based includes a small non-relativistic Schrodinger 
term characterized by a mass m and proportional to 
the square of the momentum. In addition, we include 
a dominant relativistic linear-in-momentum Dirac term 




3 


with Fermi velocity vp- It has the form^i^ 
h^k'^ 

H =-^—+ hvpikx^y - kyax), ( 1 ) 

where ax and ay are the Pauli spin matrices and k is the 
momentum relative to the F point of the surface state 
Brillouin zone. As previously discussed, representative 
values for the mass and velocity are m Ri O.lm-e (with 
nie the bare electron mass) and vp of order lO^m/s, 
respectivelyi^. When a magnetic field B is oriented per¬ 
pendicular to the 2D surface states, the magnetic energy 
scales for the relativistic (Ei) and non-relativistic {Eq) 
pieces of the Hamiltonian are of order 10 meV and 1 meV, 
respectively, for B — IT. By definition, Ei = hvp ^JeB/h 
and Eq = heB/m both involve the magnetic length 
Ib = \/k/{eB). In the Landau gauge [with magnetic 
vector potential A = (0, Bx, 0)], the Hamiltonian can be 
written as 


Pure Relativistic Topological Insulator 
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H = 


mPg L 


a^a ■ 


^ t 
-nVF-j—CL 
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—hvp-—a 
Ib 


ml^ 


a'a + 


1 


where 


FIG. 1. (Color online) (a) Perfect Dirac surface states in 
(2) the pure relativistic limit, (b) Hourglass shaped spectrum 
when a subdominant non-relativistic term is present. The 
solid dots with labels 0,1^,2^,... show the energy of the LLs 
in a magnetic field. The shading indicates the conduction 
(red) and valence (blue) bands. 
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are the creation and annihilation operators of the har¬ 
monic oscillator Hamiltonian. The energy of the LLs is 


where a|A^) = y/N\N — 1) and |A1) = 

^/N + 1 |A1 -I- 1). The spin-up amplitude is 
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where s = ± for the conduction and valence band, re¬ 
spectively. A schematic comparison between the LLs of 
the pure relativistic limit {Eq = 0) and TI is shown in 
Fig. [TJ The LLs (dots) are overlaid on the respective 
B = 0 band structures. In the TI, the N = 0 level is 
no longer pinned at zero energy but has shifted upwards. 
The symmetry between the conduction (red) and valence 
(blue) levels is also lost with finite Eq. In the following 
numerical work, we will take Ei/'/B = 10. 4 meV/Tr 
and Eq/B = 1.1 meV/T as characteristic values of a TI. 
In terms of the spin-up and spin-down amplitudes, the 
eigenfunctions of the state |A^s) are 


\Ns) 


CijN) 


( 6 ) 


gives the spin-down component. 

Using the Kubo formula, the dynamical conductivity 
aaf}{£l) is given by 


1/3 (H) = 


ih 
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Jm.s' — fN,s 


B N,M=0 

S,S^ = ± 


Ns\ja\Ms'j (^Ms'ljplNs 
h£l £m,s' — £n,s T ik/ {2 ,t^ 


( 9 ) 


where fpr.s is the Fermi function for state N in band 
s, a and j3 index the x and y components, t is a phe¬ 
nomenological optical scattering time taken to be a large 
constant and ja = eVa = {e/k){dH/dk^) is the current 
operator. At zero temperature, //v^s can be replaced by 
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the Heaviside step function 0(/z — £n,s)'-, here, ^ is the 
chemical potential. The necessary velocity operators are 


fi , » . h — a ^ 

Vx = —kx + VFCfy = l -^=-- \-VF(Jy, (10) 


m 


and 


Vy = 


eBx 


- VfCTx = 


m ^/21 b 

h o' -\- a 
\/21b 


-vf^x- (11) 


Calculating the appropriate matrix elements, we obtain 
the real part of the longitudinal conductivity^! 



X _ H _ 

(H + £m,s' ~ £n,s)^ + 

X [J-{Ns;Ms')dN,M-i + Ns)dM,N- 


where rj = h/ (2r) is a constant phenomenological optical 
scattering rate and!! 


T[Ns-Ma') = 


X 


n'\ ni 




Eq 

” V2Ei 

N,s + VlV + 1C 


4" ^ 

M,s'^ 



(13) 

n 2 


The imaginary part of the conductivity is found by re¬ 
placing the r] factor in the numerator with fl -|- £m,s' — 
£n,s- Similarly, the real part of the transverse Hall con¬ 
ductivity is!! 


Re 


\ ^xy(S^) 
\ e^/h 



CXD 

E 

iV,M=0 

s,s'=± 


fM,s' — fN,s 
£n,s £m^s' 


(14) 



H - l - £m,s' ~ £n,s 

(H -l- £m,s' ~ £n,s)‘^ + 

X [J^(A^s; Ms')Sn,m-i — E{Ms'; Ns)5m,n-i] 

Here, the imaginary part is given by the replacement H-|- 
£m,s' — £n,s —>■ —V in the numerator. 


FIG. 2. (Color online) Real part of the longitudinal optical 
conductivity as a function of photon energy at charge neu¬ 
trality (/i = 0) for (a) the pure relativistic case and (b) with 
the addition of a non-relativistic contribution. Three values 
of B are considered: B = O.IT, B = IT and B = 5T. 


III. NUMERICAL RESULTS FOR THE 
DYNAMICAL CONDUCTIVITY 

To set the stage for understanding the semiclassical 
limit, we start by examining how the dynamical conduc¬ 
tivity evolves as the magnitude of the magnetic field B is 
progressively decreased towards B = 0. As we are inter¬ 
ested in modifications to the semiclassical limit brought 
about by the subdominant non-relativistic (Schrddinger) 
contribution in the Hamiltonian [Eqn. ([IJ], differences 
between the pure relativistic (Dirac) and TI are empha¬ 
sized. 

In Fig. m we compare Re(Ta;a;(D) for a typical TI with 
the ideal relativistic system. The top frame is for Eq = 0 


(i.e. no Schrodinger mass term). In both cases, we take 
the Dirac magnetic energy scale to be Ei/'/B = 10.4 
meV/-\/T and set the chemical potential at charge neu¬ 
trality (/r = 0). Three values of magnetic field are 
considered: B=0.1T (solid black curve), B=1T [dashed 
green/blue curve for (a)/(b)] and B=5T (dash-dotted red 
curve). Therefore, Ei = 3.29 meV, 10.5 meV and 23.3 
meV, respectively. It is these energies that determine 
the energies of the various LLs. For Eq = 0, Eqn. dSD 
gives = s\/2NEi for N = 1,2,3,... and £o = 0 
for V = 0. In the pure relativistic case, one has per¬ 
fect particle-hole symmetry about the neutrality point. 
In the numerical work, we have taken the optical scat¬ 
tering rate rj of Eqn. m to be 0.5 meV which provides 
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a width to the optical absorption lines that otherwise 
correspond to Dirac (5-functions. First, we note that a 
universal interband background becomes well defined at 
high photon frequencies in the solid black curve. In units 
of this takes the value of 1/16. For larger mag¬ 

netic fields (dashed blue and dash-dotted red curves), the 
energy spacing between the optical peaks has increased 
and the background is no longer resolved. In all cases, 
the first line is an intraband transition from = 0 to 
(iV, s) = (1, -I-) while all other lines are due to interband 
transitions from {N, —) —>■ (TV ±1, -I-) with iV = 1, 2, 3,.... 
The two interband lines for a given N are degenerate in 
energy for a given value of N. These results are for com¬ 
parison with those presented in the lower frame which 
apply to the case when the Schrddinger energy Eq is not 
set at zero {Eq/B = 1.1 meV/T). For B = O.IT, Eq is 
small compared to Ei. For B = 5T, it has increased to 
5.5 meV which is roughly a factor of four smaller than 
El. Still, the energies at which the various lines occur 
have not shifted a great deal when compared to the top 
frame; however, we can clearly see that the intraband line 
for B = 5T has shifted downward upon the introduction 
of a Schrodinger term. As is clearly seen in the dashed 
blue curve, for B = IT, the interband transition peak has 
split into two with the energy of the splitting given by the 
value of Eq. This can be traced to the fact that, in a TI, 
the lines {N, ^ {N + 1, -I-) and {N, ^ {N — 1, -I-) 

are no longer degenerate in energy. As B is reduced in 
magnitude, the splitting of the interband peaks into dou¬ 
blets is no longer resolved (see the solid black curve for 
B = O.IT). This is because Eq itself becomes small and 
the residual scattering rate that we have used (ry = 0.5 
meV) washes out the splitting. The universal background 
seen at higher values of photon energy D rapidly becomes 
constant and equal to e^/(16?i), completely independent 
of the Schrddinger mass term. In fact, because we have 
taken yt = 0 (charge neutrality), the entire curve becomes 
constant independent of 11 for i? = 0. 

It is illuminating to also consider the evolution of the 
Hall conductivity as H —>■ 0. We know that at i? = 0, 
this quantity [Re(T 2 ;y(f 2 )] must be zero rather than the 
constant value seen above for Re(T 3 ,a,(ri). In Fig. [31 we 
provide a comparison between the pure Dirac and TI for 
the real part of the Hall conductivity. The same three 
values of magnetic field are considered and the notation 
is the same as for Fig. [3] Again, the energy of the first 
peak/valley is only slightly affected by the introduction 
of an Eq term. However, there are important differences 
which must be noted. Namely, in the top frame, the 
peak is always first followed by a valley at slightly higher 
energy. In the bottom frame, it is the opposite. Asso¬ 
ciated with this reversal is the fact that the DC limit 
of ReCT 2 ,y(H) which gives the quantized Hall conductivity 
has the same magnitude (l/(47r) in our units) but is of 
opposite sign. This change has been brought about by 
the introduction of the Schrddinger term Eq; this does 
not change the quantization in the plateau height (which 
remains relativistic) but does change the energy (chemi- 



FIG. 3. (Color online) Real part of the transverse Hall con¬ 
ductivity for the same parameters as Fig. O 

cal potential) at which the step to a new plateau occurs. 
For a TI, the step from —1/2 to 1/2 quantization (in 
units of e^/h) occurs at finite /r rather than at n — 0^. 
This is due to the first LL sitting at an energy Eq/2 
rather than at zero. Some further explanation of our 
language is required. As we have stated, in our units, 
the universal interband background has a value of 1/16 
or Re(Ta;x(II) = e^/{16h) = e^7r/(8h) which is a factor of 
four smaller than the well known value of e'^TTl(2h) for 
graphene. In graphene, this value accounts for a valley- 
spin degeneracy which is not present in a TI. Further, 
the Hall plateau quantization is ±1/2, ±3/2, ±5/2,... in 
graphene (in units of e^/h) if the four fold degeneracy is 
removed. Another striking difference between the pure 
relativistic [frame (a)] and the TI case [frame (b)] is the 
additional LL structure seen in the lower frame which is 
absent in the upper frame. This is traced to the split- 
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ting of the optical transitions between N — 1 ^ N and 
N —>■ A^ +1 for a TI. These oscillations are clearly signifi¬ 
cant at larger valnes of B but are progressively reduced as 
B becomes small. Note that for the B = O.IT (solid black 
curve), RetTa;y(n) ^ 0 for most of the il range shown. At 
B — 0, the entire curve is of course zero. 

So far, we have considered /i = 0 to get a first under¬ 
standing of the differences and unchanged characteristics 
of the magneto-optics which are brought about by the 
introduction of a subdominant Schrbdinger piece in the 
Hamiltonian. Next, we will concentrate on finite fj,. In 
fact, our interest will be on the semiclassical regime for 
/i much larger than the magnetic scales Ei and Eq. In 
Fig. m we show the results for the real parts of the longi¬ 
tudinal conductivity Recr 2 ,a;(r 2 ) vs. H [frame (a)] and the 
Hall conductivity Recra;i/(f^) vs. H [frame (b)]. The mag¬ 
netic field is set at B = 0.2T and three values of /r are 
considered. The solid black curve corresponds to /r = 0 
and is for comparison. The dashed red curve applies to 

= 8 meV and the dash-dotted green line is for ^ = 10 
meV. We want to focus on the first intraband peak and 
understand its evolution as n is increased. It is this line 
which eventually evolves into the cyclotron resonance line 
of the semiclassical theory. We note that, as p. increases, 
the position of this peak moves to lower energies and ac¬ 
quires added intensity (optical spectral weight). Frame 
(b) gives ReCT 2 ,y(H). We note that, in these curves, the 
position in energy of the first valley-to-peak (for ^ = 0) 
and peak-to-valley (for finite /r) structures also evolves 
towards lower energies and tracks the intraband peaks in 
Reo’xa:(f^) of the top frame. In particular, we emphasize 
the DC value of Recra;y(D) which gives the Hall plateaus. 
For fj, = 0, it is —1/2 [i.e. —e^/(47rS)] while for ^ = 8 
meV it has moved to 3/2 and then to 5/2 (in units of 
e^/h) for fj, = 10 meV. These values are independent of 
the mass term in the Hamiltonian. We now turn to the 
semiclassical regime fi Ei and Eq. 



FIG. 4. (Color online) (a) Real part of the longitudinal AC 
conductivity as a function of photon energy for B = 0.2T and 
three values of chemical potential: fj, = 0 (charge neutral¬ 
ity), /i = 8 meV and p = 10 meV. (b) Real part of the Hall 
conductivity for the same parameters as (a). 


IV. SEMICLASSICAL LIMIT 


For definiteness, we assume p > 0 (i.e. the conduction 
band is partially filled). In this case 


We begin with Eqn. (jT^ for the real part of the AC 
longitudinal conductivity. We work in the limit of p 
much larger than the magnetic quantization such that the 
chemical potential falls between the and {Nc + 1)*^ 
LLs with Nc very large and given by 


f CTxxi^) \ _ ^ fN,+ — fN-l,+ 

\ e^jh J 27r £]v-i,+ - £n,+ 

yE(N-l+;N+) 

(fi -l- £n,+ — £n-i,+)^ + 
7jE(N-1+;N+) ] 

(D -I- £jv-i,+ — £n,+)^ +‘> 7 ^ J 



(16) 


p — EqNc + 


l2NcEf 



From here on, the -I- index shall be implied (and sup- 
(15) pressed) as we are only interested in s,s' = -I-. In our 
limit, the first factor of the sum becomes —df^/dSN and 
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the energy difference £n—£n-i = WcR(/r) is the cyclotron 
resonance energy for /x > 0. At zero temperature, 


Re 


' ■' N =0 


For B 


(n + Wcr)^ +77^ {£ 1 - Wcr)^ + 

0, Aq <C Ai and Eqn. (113 gives 


Eq 
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2/i 


- 1 


(18) 


The sum over N in Eqn. m also changes into an inte¬ 
gral. Therefore, the Dirac ^-function S(fi — fjv) means 
that the entire integrand is to be evaluated at iV = iVc. 
This gives 
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( 21 ) 
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(D-I-Wcr)^-I-(D - Wcr)^-I- 


The only £l dependence that remains is in the curly 
bracket and the integral over photon energy of this term 
is normalized to a value of one. It provides an analytic 
formula for the frequency profile of the cyclotron reso¬ 
nance line. The remaining factors in Eqn. dan give the 
optical spectral weight of the cyclotron resonance which 
we denote Wcr, where 


WcR = 




2 m _ Y 


h 2 


1-k 


2m 


-E{Na-l\Nc). (22) 


We now need to evaluate F{Nc — 1; W)- Since is very 
large as R —>■ 0, we can replace A^c — 1 by to get 


^(iVc; N,) 


M pi _ ^0 


(23) 




which follows from Eqn. (nsi). Using Eqns. © and ®, 
we note that 


and 


'^Nc 


pi 


Eol2 

2y/W~^Ei 


(1 , Eo/2 


I 2 2y/WcEi ■ 

Using Eqn. (ITOl) for \/2Nc, we find 

Eq Eq 1 


^ 2 NcEi mvj, h I 2m _ 2 

V 


(24) 


(25) 


(26) 


which goes to zero as R —>■ 0. Therefore, ~ —\f\j2 
and ~ Thus 


1 Eo 1 




F{Np,N,) 


When this is substituted into Eqn. (I22L we obtain 


WcR = 


e mVf 

h 8 


2/i 


- 1 


2/i 


(27) 


(28) 


for the optical spectral weight of the cyclotron resonance 
line. Including negative chemical potentials, we find 


WcR = 


h 8 


2/i 


mv7 


- 1 


2/i 


mv 


2 ’ 


(29) 


for /i ^ 0. The profile in photon energy of the line is 
given by 


77/7 


(D -I- Wcr)^ -I- 7/2 (D - Wcr)^ + V 


rjl? 


(30) 


for both /i ^ 0. This formula is unchanged from the result 
found for graphene^ except that the cyclotron frequency 
Wcr has been altered. By definition, the cyclotron res¬ 
onance energy Wcr at chemical potential /t is given by 
~ ^Wc-i = <^cr in the limit of large N^. determined 
by /i. We get 



, (31) 


, , heBvp (^ I 2/i 


which does depend on the Schrodinger mass. In Eig.[5])a), 
we plot the optical spectral weight (Wcr) of the cy¬ 
clotron resonance line in a TI as a function of |/i|. Three 
cases are considered. The dash-dotted green curve is for 
comparison and corresponds to the pure relativistic limit. 
In this instance, Wqr increases linearly with |/i|. For a 
TI with finite Schrodinger contribution, the red curve 




























































FIG. 5. (Color online) (a) Optical spectral weight under the 
cyclotron resonance line as a function of |/r| for the pure rela¬ 
tivistic and TI regimes, (b) The corresponding results for the 
cyclotron frequency at B = O.IT. 


(/i > 0) is always above the relativistic case while the 
dashed blue curve < 0) is always below it. In frame 
(b), we present the complimentary results for the |/i| de¬ 
pendence of the cyclotron resonance frequency Wcr for 
B = O.IT. Again, we see deviations from the pure rela¬ 
tivistic limit with the solid red curve (/i > 0) and dash- 
dotted blue curve (/r < 0) falling above and below this 
case (dash-dotted green), respectively. It is interesting to 
obtain a first correction to Wcr and Wcr as a function of 
1 /to. Trivially, the to —>■ oo limit of Eqn. (1311) gives the 
well known result for Wcr in the pure relativistic regime; 


that is, Wcr = heBvp/\n\. The first correction is a factor 
of 

which, as seen in the Fig. [S](b), moves the curve up for 
fi > 0 and down for ^ < 0. In Fig. [Ub), only |^| > 10 
meV is shown as cyclotron resonance is defined for large 
/r. For the optical spectral weight, the first correction 
gives 

Again, the additional contribution to the relativistic limit 
is positive (negative) for ^ > 0 (/x < 0). 

So far, the results have been presented as a function of 
chemical potential. Of course, the electronic density of 
states in a TI is not symmetric with respect to the neu¬ 
trality (Dirac) point. It rises above the pure Dirac limit 
for negative energies and is reduced below this regime for 
positive fj,. This is shown in the inset of Fig. [51 Thus, for 



FIG. 6. (Golor online) The number of dopants (n) as a func¬ 
tion of the magnitude of the chemical potential. The pure 
relativistic case is contrasted with that when a subdominant 
Schrodinger term is included. Inset: density of states as a 
function of oj. The shading identifies the conduction (red) 
and valence (blue) bands. 


a given magnitude of /x, there are more holes {n-h) than 
electrons (ue) involved. It is of interest to see how large 
the asymmetry in Wcr and Wcr are when plotted with 
respect to an equal value of doping (holes or electrons). 
To accomplish this conversion, we use the relationship 
between /x and n. Namely, 




( 34 ) 
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Using Eqn. (m, we find 


WcR = e' 


, \/47m I 2 ^ \/47rn \ 

16 y 2mVp J 


(35) 


and 


WcR — 


eBhvp 


2mv'% j 


\/47rn \ 
2mv'^ j 


(36) 


for the optical spectral weight and frequency of the cy¬ 
clotron resonance line, respectively. In both cases, ± 
corresponds to electrons and holes, respectively. Here, 


l/^l 


mvp 

2 


•\/47rn 

2 

^ ^ 2mv'jp 

-1 


(37) 


The main frame of Fig. [5] shows the /i dependence of 
the density n in units of hr'^Vp^ meV^ for the pure rela¬ 
tivistic case (dash-dotted green) and TI with ^ > 0 (solid 
red) and /r < 0 (dashed blue). The intersection of the 
curves with the dotted vertical line at |/r| = 20 meV (cir¬ 
cles) emphasizes the different electron densities at fixed 
fj,. The horizontal dotted line gives the different chemical 
potentials (marked by a x) associated with a fixed parti¬ 
cle number (n = 200 in our units). As expected, for fixed 
/i, there are more holes than electrons. To get a constant 
number of dopants, the magnitude of the chemical poten¬ 
tial needs to be considerably larger in the electron case 
than for holes. The impact of this conversion from chem¬ 
ical potential to dopant density is shown in Fig. [7] for 
the optical spectral weight [frame (a)] and cyclotron res¬ 
onance frequency [frame (b)]. Large differences between 
electrons and holes are noted. 

Next, we consider how the optical spectral weight un¬ 
der the intraband line (lUintra) [i.e. the area under the 
first peak of Fig. |4l[a)] and its position in energy (wintra) 
at finite B approach the cyclotron resonance results as B 
is decreased. This is shown in Fig. [8l In the top frame, 
we reproduce as the dotted red curve, the results previ¬ 
ously obtained in Fig. [5] (solid red) for the spectral weight 
of the cyclotron resonance line Wcr in units of e^/fi meV 
as a function of chemical potential for /i > 0. This curve 
is to be compared with the related results for the opti¬ 
cal spectral weight Wintra under the intraband line. The 
dashed blue curve is for B — 0.5T while the solid black 
curve corresponds to B = O.IT. As the magnetic field 
is reduced and p, is increased, our results for Wintra ap¬ 
proach the dotted red curve for Wcr as we expect. The 
two curves should agree in the limit of /r much larger 
than the magnetic energy. For B — O.IT, this condition 
is /r » 3.3 meV. In Fig. [51[b), we show the energy of 
the intraband transitions Wintra in meV as a function of 
/J, for the same two values of magnetic field: B = 0.5T 
(dashed blue) and B = O.IT (solid black). The dotted 
green and red curves are for comparison and give Wcr vs. 
/i for H = 0.5T and O.IT, respectively. As the value of /i 



FIG. 7. (Color online) (a) Optical spectral weight under 
the cyclotron resonance line as a function of the number of 
dopants n for the pure relativistic and TI regimes. Inset: den¬ 
sity of states as a function of uj. The valence and conduction 
portions are given by the blue and red shading, respectively, 
(b) The corresponding results for the cyclotron frequency at 
B = O.IT. 


is increased, the condition that /i be much greater than 
the magnetic energy scale is better satisfied and Wintra 
merges with Wcr. 


V. THE INTERBAND BACKGROUND 


In the previous section, we emphasized the intraband 
contribution to the optical conductivity. However, the 
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FIG. 8. (Color online) (a) Optical spectral weight under the 
intraband line as a function of /r for a TI. Two values of 
magnetic field are shown: B = O.IT and B = 0.5T. The 
dotted curves give the corresponding results for the cyclotron 
resonance spectral weight, (b) The position in energy of the 
intraband line for the same parameters as above. 


interband part is also of interest and corresponds to all 
but the first magneto-optical lines of Figs. [UHl At larger 
values of photon energy aird in the small B limit, these 
provide the universal background which is given by the 
known formula 

Re(Ta:a;(f^) = ~ nrax[2A, WmiJ), (38) 

where Wmin is the threshold energy for the onset of inter¬ 
band transitions. In Ref4i, this formula is derived from 
the form of the conductivity in the case of i? = 0 which 
involves an integration over momentum. Eqn. (1381) prop¬ 


erly reduces to the graphene resul t which differs from 
this equation only in that Wmin = 2/r. In the appendix, 
we provide an alternative derivation of this formula. We 
start with Eqn. m for the longitudinal conductivity at 
finite B and take the limit of i? —^ 0. This serves as 
a check on our work. The inset of Fig. IH] identifies the 
energy Wmin for a particular value of |/i|. It is smaller 



FIG. 9. (Color online) Energy cutoff of interband transitions 
as a function of |/i|. Positive and negative are considered. 
The results are compared to the pure Dirac limit. Inset: 
sketch of the dispersion curves. The arrows identify the inter¬ 
band transitions that define Wmin. Blue corresponds to /r < 0 
while red is for > 0. 


for ^ > 0 (red arrow) than for /j, < 0 (blue arrow). The 
value of Wniin is determined by £+{kc) — £-{kc) with kc 
corresponding to p, = £+{kc) for /i > 0 and /i = £-{kc) 
for ^ < 0. Here, £±{k) = h^k"^ /{2m) ± hvpk. This gives 


Wmin(Ai) = 2.mv% 



(39) 


which reduces to the known value of 2/r. in the pure 
relativistic limit. For a TI with a subdominant non- 
relativistic contribution, our results for tLimin vs. |^| are 
presented in the main frame of Fig. |9l The dash-dotted 
green curve is for comparison (m —>■ oo). For /i > 0, the 
curve (solid red) falls below the green line. For /r < 0, 
Wmin (dashed blue) is above the relativistic limit. As 
we expect, all three curves come together at small |/r| 
since, in that limit, the relativistic linear-in-momentum 
magnetic energy dominates the non-relativistic quadratic 
term. As |^| is increased, deviations from the pure rela¬ 
tivistic result become significant. 

In Fig. m a), we show both the intraband line (Drude- 
like peak near H = 0) and the constant interband back¬ 
ground which sets in at il = Wmin- The absolute value of 
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FIG. 10. (Color online) Real part of the longitudinal AC 
conductivity as a function of photon energy, (a) \fj,\ set at 
30 meV. A different interband onset is seen for /r ^ 0 when 
the Schrodinger term is finite, (b) Corresponding results for 
a fixed dopant concentration nh^Vp = 200 meV^. 


^ is set at 30 meV so that the various onsets for ^ > 0 
(solid red) and /x < 0 (dashed blue) in a TI, and the pure 
Dirac curve (dash-dotted green) are clearly separated. 
The height of the background remains at (l/16)e^//i in¬ 
dependent of the mass term. In frame (b), we show re¬ 
sults for fixed dopant density n set at 2QQh~'^Vp^ meV^. 
Again, the onset of interband transitions is changed by 
the subdominant Schrodinger term. Note that these re¬ 
sults were obtained from finite magnetic field calculations 
for B = O.IT and consequently small wiggles from the 
underlying LL structure remain. These wiggles are most 
pronounced for small values of photon energy. Since we 
have included a residual scattering rate rj = 0.5 meV, the 
Drude term centred at D = 0 is broadened. It would be a 


Dirac 5-function in the clean limit. This smearing param¬ 
eter is also responsible for the small but finite conductiv¬ 
ity seen between the Drude and interband contributions. 
Increasing the residual scattering rate would increase the 
value of this finite conductivity and reduce the small wig¬ 
gles seen in the interband transitions of Fig. [TUI 


VI. DISCUSSION AND CONCLUSIONS 

We start with the Kubo formula for the real part of 
the longitudinal and transverse (Hall) conductivities in a 
magnetic field B which is oriented perpendicular to the 
surface of a TI. The magnetic field creates LLs which lead 
to absorption peaks in the real part of the longitudinal 
optical conductivity. The introduction of a subdominant 
Schrodinger contribution to the Hamiltonian is known to 
split each interband peak into a doublet. Here, we find 
that there are associated peaks in the real part of the 
dynamical Hall conductivity ReCTa:y(D) which are entirely 
due to the non-relativistic term. 

A primary interest is the semiclassical regime which 
emerges when the magnitude of the chemical potential 
|/x| is much larger than both the relativistic (i?i) and non- 
relativistic {Eq) magnetic energy scales. Nevertheless, we 
have found it illuminating to begin with a study of the 
evolution of both the dynamical conductivities Re(Ta;a;(D) 
and Recr 2 ;y(D) as the magnitude of the magnetic held 
is gradually reduced towards the zero-held limit. For 
B — 0, the intraband optical transitions at hnite doping 
provide a Drude contribution to Re(Ja;a;(D) and the inter¬ 
band transitions give an additional universal background 
the height of which is independent of the Schrodinger 
mass; however, its onset is not. This step is also not 
particle-hole symmetric whether it be as function of /x or 
the number of dopants (n). In the process of consider¬ 
ing the interband background, we have provided a new 
derivation of the relevant formula. The derivation pro¬ 
ceeds from our formula for Recra;a;(D) at hnite magnetic 
held and the limit of H —>■ 0 is formally applied. At hnite 
H, we hnd LL structures associated with the interband 
transitions in Recr 2 :y(D) which exist only when m is hnite. 
They are gradually reduced in amplitude with decreasing 
B; for B = 0, the entire Hall conductivity is zero. For 
hnite doping, the intraband transitions provide a peak 
in Re(Ta;a:(fl) at low photon energy D which evolves into 
the cyclotron resonance line as B gets small or |/x| gets 
large (such that |/i| is much greater than Ei and Eq). We 
hnd a corresponding peak-valley structure in Re(Ta;y(D) 
with its DC limit equal to the Hall plateau quantiza¬ 
tion which keeps the pure relativistic limit sequence of 
(±1/2, ±3/2, ±5/2, ...)e^jh. The value of /x at which the 
transition from one plateau to the next occurs is affected 
by the non-relativistic mass term. One consequence of 
this fact is that at charge neutrality (/x = 0"^), the DC 
Hall conductivity (in units of e^/h) becomes —1/2 for 
Eq 0 while it is 1/2 for Eq = 0. This also implies that 
the structure seen in Re(T 2 ;j,(D) at small D is switched 
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from valley first to peak first. 


When the chemical potential is large compared to the 
Dirac and Schrodinger magnetic scales, the intraband line 
moves toward the semiclassical cyclotron resonance line. 
We provide a simple analytic formula for the cyclotron 
resonance frequency, its optical spectral weight and its 
line shape which is found to be unaltered from the pure 
relativistic case. However, both the cyclotron frequency 
(wcr) and optical spectral weight {Wcr) are changed. In 
addition, due to the Schrodinger-mass induced particle- 
hole asymmetry, both Wcr and Wcr are different for 
positive and negative fi. Differences are also evident 
for a fixed doping of electrons or holes. We provide a 
new derivation of the universal interband background of 
the longitudinal conductivity. The work is based on the 
Kubo formula at finite magnetic field and a formal limit 
of H —>■ 0. We recover the result obtained when no mag¬ 
netic field is considered. The height of the background is 
completely independent of the Schrodinger mass m but 
its onset at fl = Wmin is not. This onset not only depends 
on m but also depends on the sign of the chemical po¬ 
tential or (for a fixed value of dopant density) whether 
they are electrons or holes. 


Appendix 


In this appendix, we provide a derivation of the uni¬ 
versal background formula [Eqn. (jSS])] as the H —?> 0 limit 
of Eqn. (HH). The two optical matrix elements of interest 
are for transitions from {N — I, —) to (iV, -I-) and {N, —) 
to (iV-l-1, -|-). In addition, these are only needed for large 
N with N = Nc given by /i = ^N„+- For simplicity, we 
have assumed /r > 0. The negative ^ case follows from 
similar considerations. We can take 77 —>■ 0 in Eqn. (HU) 


to get 



(A.l) 


I fN+l,+ + 1 + )5(D + - £n,-) 

I ON,- — CAr+l,+ 

+ /jV-l,+ -/jV,- ^(^_; ^ ^ J ^ 

ON,- — ON-l,+ J 


In the limit 5—^0 and large fi, the Dirac ^-functions 
in Eqn. (|A.1I) restrict N to values near Nc where H = 
2\/2NcEi which also needs to be larger that Wmin- The F 
factors can be replaced by F{Nc—; Nc+) and the energy 
denominators by —2\/2NcEi since Eq goes to zero as 
H ^ 0. Finally, the thermal factor with a — needs to be 
replaced by 1 and the one with -|- by 0. This allows only 
transitions with an energy larger than Wmin ot N > Nc- 
Thus, we have 



I 

V^cEi 


F{Nc-;Nc+) 


X 6{£l + £n,+— Sn,-), (A.2) 


with H > Wniin- The sum over N can be replaced by an 
integral over a continuous variable N. Integration is then 
preformed on the Dirac 5-function which gives a factor 
oi [{d/dN)2'/2NEi]~^\N,,- Therefore, 


2 £;i 

(A.3) 

But, F{Nc-; Nc+) = 1/4; so, 

Re(Ja;a;(H) = — (A.4) 


Re 


CTxxj^) \ _Ei 


VS i 2 y/m^E: 


-F{Nc-;Nc+) 


which is completely independent of the Schrodinger mass 
TO. However, this constant background starts at = 
Winin and so, as we have seen, this onset does depend on 
TO as emphasized in Eqn. (1551) . 
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